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Residual correlations between decay products due to a combination of both cor- 
relations between parents at small relative velocities and small decay momenta are 
discussed. Residual correlations between photons from pion decays are considered as 
a new possible source of information on direct photon fraction. Residual correlations 



in p7 and pA systems due to pT> interaction in final state are predicted based on 
the pTP low energy scattering parameters deduced from the spin-flavour SUg model 
by Fujiwara et al. including effective meson exchange potentials and explicit flavour 
symmetry breaking to reproduce the properties of the two-nucleon system and the 
low-energy hyperon-nucleon cross section data. The P7so residual correlation is con- 
centrated at k* ~ 70 Mev/c and its shape and intensity appears to be sensitive to the 
scattering parameters and space-time dimensions of the source. The pA s o residual 
correlation recovers the negative parent pTP correlation for k* > 70 Mev/c. The 
neglect of this negative residual correlation would lead to the underestimation of the 
parent pA correlation effect and to an overestimation of the source size. 



I. INTRODUCTION 



Due to the effects of quantum statistics (QS) and final state interaction (FSI), the mo- 
mentum correlations of two or more particles at small relative velocities, i.e. at small relative 
momenta in their center-of-mass (cm.) system, are sensitive to the space-time characteristics 
of the production processes on a level of fm = 10~ 15 m. Consequently, these correlations are 
widely used as a correlation femtoscopy tool providing a unique information on the reaction 
mechanism which is hardly accessible by other means (see, e.g., recent reviews l|, |2j). 

The momentum correlations of two particles with four-momenta p\ and P2 are studied 
with the help of the correlation function H(pi,P2) which is usually defined as the ratio of 
the measured distribution of the three-momenta of the two particles to the reference one 
obtained by mixing particles from different events of a given class, normalized to unity at 
sufficiently large relative momenta. It can be also written as the ratio of the two-particle 
production cross section to the product of the single-particle ones 

^ Pi, P2 = N — , 1 

era /crpi ■ a 6 a / crp2 

where N is the normalization factor which is sometimes taken weakly dependent on the 
relative momentum to account for the effect of possible non-femtoscopic correlations. 

At a first glance, one can hardly expect any correlations between photons from the decays 
of different neutral pions. The spatial separation between such photons is of the order of 
10 6 fm and the corresponding Bose-Einstein enhancement is extremely narrow and practi- 
cally unobservable. Nevertheless, due to femtoscopic QS correlations between parent pions 
at small relative momenta as well as due to a small decay momentum, correlations be- 
tween decay photons from different neutral pions should exist and have been experimentally 
observed [3|, |4j. The small decay momentum guarantees that a small relative momentum 
between photons corresponds to a small relative momentum between parent pions. For this 
kinematic reason the QS correlation between neutral pions is transferred to decay photons, 
being however smeared due to randomly distributed directions of the decay three-momenta 
in the respective parent rest frames. We shall refer to such correlations as the residual ones 

inn 

(see also [3|, |4], 15]). 

The residual correlations are important not only for two-photon system. For example, the 
two-baryon correlations are also affected by residual correlations arising from the FSI and 
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QS correlations in the systems involving strange baryons since the decay momenta in their 
decays (e.g., A — ► pir~ or S° — > A7) are not so large to destroy the original correlations jlj]. 
The relative importance of the residual correlations is however quite different for two-photon 
and two-baryon systems. Usually, the residual two-photon correlation almost completely 
dominates over the correlation of direct photons. The residual correlations in two-baryon 
and other systems become more and more important with the increasing collision energy 
due to increasing fraction of the produced strange particles. Their analysis is therefore an 
up-to-date task. 

The residual correlations do not represent only the distorting effect which introduces 
additional systematic errors in correlation studies. We would like to pay attention to the 
fact (to our knowledge, for the first time) that the residual correlation itself is a valuable 
source of femtoscopic information. 



II. TWO-PHOTON CORRELATIONS 

To study the residual correlations between photons (7^0) from neutral pion decays 
(7T° — ► 77), we have assumed that the photons are produced either through these decays or 
"directly" (similar to the production of 7r°'s) and generated neutral pions and direct photons 
(7d) according to the same thermal-like momentum distribution 

dN/dpoc(p 2 /E)exp(-E/T ), (2) 

where T =168 MeV. Usually, there is the experimental threshold in the energy of detected 
photons. To take it into account we have rejected photons with E 7 < 80 MeV. Such a cut 
increases the relative strength of the correlation of direct photons and modifies the relative 
strength and shape of residual correlations. 

The fraction d = N(jB)/[N( i j ir o)+N(ji))] of direct photons was considered as a simulation 
parameter (d = 0, 5, 10, 20%). Only the QS correlation has been introduced by giving each 
simulated pair of neutral pions and direct (unpolarized) photons a weight (see, e.g., [7() 

ft(pi,p 2 ) = 1 + A(cos(2k*r*)) 

= 1 + Aexp(-Q 2 r 2 ), (3) 

where Q = 2k* and r* are respectively the relative three-momentum of the two particles and 



spatial separation vector of their emitters in the two-particle rest frame; Q = \J —{pi — p-i)' 



for the considered equal-mass particles. The correlation strength parameter A = 1 for pions 
and 1/2 for unpolarized photons. A Gaussian redistribution with the same radius parameter 
ro has been assumed in the averaging in Eq. ([3]) over the spatial separation r* for both pion 
and photon emitters: 



This is a reasonable assumption for pions emitted with moderate transverse momenta. How- 
ever, for photons, it leads to the non-realistic dependence of the correlation function on the 
outward component of the relative momentum in the longitudinally co-moving system (the 
component in the direction of the pair transverse momentum) due to the diverging Lorentz 
factor of the transformation to the two-photon rest frame at Q —>■ 0. 

The two photon correlation functions corresponding to the radius parameter ro = 5 
fm and different direct photon fractions d are shown in Fig. [1] as functions of the relative 
momentum Q = 2k* in the two-photon cm. system assuming the ideal three-momentum 
resolution. The correlation functions are normalized to unity at Q > m w . The peak at 
Q = m^o is related to photon pairs from the same tt°. The width of the peak in a real 
experiment depends on the three-momentum resolution. The residual correlations between 
decay photons from different 7r°'s (~ 81% contribution for d = 0.1) result in a smooth 
structure at Q < 0.17 GeV. The uncorrelated background (~ 18% contribution for d = 0.1) 
arises from 7 7r o7rj-pairs. The pairs of direct photons (~ 1% contribution for d = 0.1) provide 
the interference enhancement with a width of 1/vq ~ 40 MeV/c. 

It is important that 

(i) the residual correlations represent a first order effect in the direct photon fraction d, 
to be compared with the second order effect of the interference correlations of direct 



(ii) the residual correlation effect appears to be wider than the interference effect for 7r°'s 
or direct photons. 

As a result, for some combinations of a large source size r and a small direct photon fraction 
d, the direct photon interference is practically unobservable while the residual correlations 
can still be used to measure d. 

It is important for the suggested method that the residual correlation function of decay 
photons (depending on the three-momentum spectrum of neutral pions, their correlations 




(4) 



photons; 



4 



1.025 



1.02 



1.015 



1.01 



1.005 



I 

I 
I 

■ ■■II 

I 
I 

I 

_ III). 
^ I 
I 

nil ft*. I ""nip 

l'lfti'iffi ■ i " 
~l 1 1 1 1 1 

^^^^^ 1 ■ 1 ■ 

"■ 1 ■ 1 
™^ "l ■ 1 ■ 

~ -■ 1 ■_ 

1 * f 

_ "i 1 H 1 ^^^^ L " 

~ ^ 1 imii-a i ■ j 

* l 1 ■ l_ ^^■^^fl 1 ■ ■ ■ ■ ■ L 

: . Mil 







0.05 



0.1 



0.15 0.2 

Q, GeV/c 



FIG. 1: The two-photon correlation functions calculated with the source size parameter ro = 5 fm 
and different direct photon fractions d = N (jd) / [N (^^a) + N(jq)]. The histograms correspond 
to d = (dashed-dotted), 0.05 (dotted), 0.10 (solid), 0.20 (dashed-triple-dotted). The thin solid 
histogram corresponds to d = 0.20 and the residual correlation switched off. 

and experimental conditions) could be predicted with sufficient accuracy. This is in contrast 
with the model-independent correlation measurement of the direct photon fractions in the 
experiment WA98 |4| exploiting the quadratic relation between the correlation strength 
parameter A and the direct photon fraction (valid in sufficiently narrow interval of the three- 
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momenta of the selected photon pairs) and nearly constant residual correlation function of 
decay photons at Q < 90 MeV/c. 



III. THE RESIDUAL CORRELATIONS IN hA AND hy SYSTEMS INDUCED 
BY THE CORRELATIONS IN hE° SYSTEM 

A. Kinematic considerations 

Another interesting example of residual correlations are the correlations in hA and h'y 
systems (where h is a hadron) induced by the hTP correlation at small relative h— S momenta 
2K in the hTP cm. system. One can express the momentum k* of the hadron h in the hA 
or h'y cm. system through the hadron three- momentum — K as 



k* = m h 



2 2/2 

z — m /m h 



1/2 



1 + m 2 jm\ + 2z 
uJ D (^h^s + K 2 ) p D K(u h + uy) 

Z = 2 1 2 C> (5) 

where p^ = 74 MeV/c is the decay momentum in the decay T° — > A7, m is the mass of the 
decay particle (a proton or a photon), cu D = [m 2 +p 2 ^) 1 ^ 2 is the corresponding decay energy, 
LOi = (m 2 + K 2 ) 1 / 2 are the energies of the particles % = h,TP in the cm. system and ( 
is the uniformly distributed cosine of the angle between the vectors po and K. 

Since the decay momentum in the S° — > A7 is rather small, the velocity of the decay-A 
is close to the parent (E°) velocity and so a substantial part of the hYP correlation at small 
relative h — TP velocities is transferred to the hA correlation at small h — A relative velocities. 
More quantitatively, it follows from Eq. (jSJ) that the parent correlation effect at K ~ of a 
width AK yields the hA residual correlation effect at (k*) ~ pD m /i/(wii + m k) of a width 
comparable with AK. Actually, for pd K m,, m, the momentum k* is practically 
independent of (: k* if(m/ms)(m/, + m^)/(mh + m.) and so, for m « ~ the 
residual correlation function recovers the parent one for pd <C fc* m. 

As for the transfer of the /zE° correlation at small relative velocities to the h'y correlation, 
the latter is shifted to (k*) = Pd/(1 + 2pd/^/i)^ 2 with the relative width Ak*/(k*) w 
AK/[\/3)U/ik(1 + PdA^/i)]) "where //^ = m h m T.I '( m h + m s) is the reduced mass of the hTP- 
system. For example, if the hadron h were a pion or a proton, the residual correlation effects 
would be respectively situated at k* w 34.8 and 68 MeV/c with the corresponding relative 
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widths AK I (332MeV / c) and AiT/(981MeV/c) much smaller than unity provided that the 
parent correlation width AK is less than ~ 100 MeV/c. 



B. Single-channel approach 

In the following we will take the hadron h to be a proton. We thus have to calculate the 
FSI correlation functions for the systems ab = pA and pTP (the FSI between direct photons 
and protons can be neglected). The two-particle correlation function at small k*- values is 
basically given by the square of the wave function of the corresponding elastic transition 
ab — > ab averaged over the distance r* of the emitters in the two-particle cm. system and 
over the particle spin projections {tJ: 



ft(Pi,P 2 ) = (IV 



= 1 + 



s 



I 2 ) 




1 


f S (k*) 


2 


r 



+ J r- K J F 1 (Qr ) - 1 ' F 2 (Qr ) 
V 7rr o r o 



(6) 



where F\(z) = Jq dxe x2 ~ z2 / z and F^^z) = (1 — e~ z2 )/z and ps is the emission probability 
of the two particles in a state with the total spin S; we assume the emission of unpolarized 
particles, i.e. po — 1/4 and p% — 3/4 for pairs of spin- 1/2 particles. The analytical expression 
in Eq. corresponds to the Gaussian redistribution (j3J). It implies a small radius of the 
FSI interaction as compared with the characteristic separation of the emitters in the two- 
particle cm. system. The non-symmetrized wave function describing the elastic transition 
can then be approximated by a superposition of the plane and spherical waves, the latter 
being dominated by the s-wave, 

V-k* ( r ) = exp(-zk r ) + / (k ) . (7) 

The s-wave scattering amplitude 

m = ^exp(2^)-l = (1/yS _. rrl; (8) 

where < rj s < 1 and 5 s are respectively the elasticity coefficient and the phase shift, 
K s is a function of the kinetic energy, i.e. an even function of k*. In the effective range 
approximation, 

1/K S = l/a s + ^d s k* 2 , (9) 
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where a s and d s are respectively the s-wave scattering length and effective radius at a given 
total spin S; in difference with the traditional definition of the two-baryon scattering length, 
we follow here the same sign convention as for meson-baryon or two-meson systems. 

One can introduce the leading correction 0(\a \ 2 d /Vq) to the correlation function in Eq. 
([H]) to account for the deviation of the wave function [j7]) from the true solution inside the 
range of the two-particle strong interaction potential 7J: 

Aft( Pl ,p 2 ) = -(V^r V£psl/ 5 (ni 2 ^(n (10) 

s 

where the function d s (k*) = 2%ld(K s )~ l /dk* 2 ; d s (0) is the effective radius. 

It should be noted that the two particles are generally produced at non-equal times in 
their cm. system and that the wave function in Eq. ([6]) should be substituted by the Bethe- 
Salpeter amplitude. The latter depends on both space (r*) and time (£*) separation of the 
emission points in the pair rest frame and at small \t*\ coincides with the wave function ip s 
up to a correction 0(\t* /mr* 2 \), where m is the mass of the lighter particle. It can be shown 
that the equal-time approximation in Eq. ([6]) is usually valid better than to few percent 
even for particles as light as pions 7, s|. 

In this paper we use the low-energy scattering parameters for hyperon-nucleon systems 
obtained within the spin-flavour SU6 quark model including effective meson exchange po- 
tentials and explicit flavour symmetry breaking of the quark Hamiltonian to reproduce the 
)roperties of the two-nucleon system and the low-energy hyperon-nucleon cross section data 

The i^-function and the low energy scattering parameters are real in the case of only 
one open channel as in the near threshold pA scattering. For pA system, we use the values 
from Table 6 of Ref. [6j: a = 2.59 fm, a 1 = 1.60 fm, d° = 2.83 fm and d 1 = 3.00 fm. 

For pTP system near threshold, there are two more open channels, nS + and pA ones, 
so, in principle, one has to solve the three-channel scattering problem. Assuming isospin 
conservation, this problem reduces to the single-channel one for isospin I = 3/2 and to 
the two-channel one for isospin I = 1/2. The NT,(I = 3/2) scattering parameters af and 
df are given in Table 6 of Ref. [6|: a^ 2 = 2.51 fm, a\j 2 = —0.73 fm, d\j 2 = 4.92 fm 
and d\j 2 = —1.22 fm. The coupling between the channel NT,(I = 1/2, S = 0) and the iVA 
channels appears to be quite weak, i.e. the elasticity coefficient rj°(k*) = 1, so the low-energy 
scattering parameters for the iVS(J = 1/2, 5* = 0) are real; in accordance with Eq. ([8]), the 



it of the energy dependence of the NTi(I = 1/2, S = 0) phase shift in figure 15 (fss2) of Ref. 
f| yields a\j 2 = —1.1 fm and d\r 2 = —1-5 fm. The situation is quite different for the channel 
NT,(I = 1/2, S = 1) which appears to be strongly coupled with the NA channels (due to 
the pion exchange potential) as demonstrated by figures 15 and 31 of Ref. fl. As a result, 
the low-energy scattering parameters for the system NT,(I = 1/2, S = 1) acquire imaginary 
parts: a\^ 2 = (—1.1 + i4.3) fm, dy 2 = (—2.2 — i2A) fm. To get these values, we have used 
the fact that the iT-function is even in k* and employed the expansions corresponding to 
the effective range approximation in Eq. (Q: 

r)(k*) = 1 - 2%ak* + 2(3a) 2 F 2 + [2(Ka) 2 Sa - 2(3a) 3 - Qd{Qa) 2 



+Qd{^a) 2 + 2MaSa 



,*3 



-4(Ka)(3a) 2 + 2(3a)' 



-4M£a(3a) 2 + 23d(3a) 3 - 2%d($la) 2 Sa 



< *4 



5(k*) = 5(0) + 3lak* + [-^(Ka) 3 + 3?a(3a) 2 + ^d(Sa) 2 

+3d9?a3a - $td($ta) 2 ] k* 3 + [-2(Ka3a ((Ka) 2 + (3a^ 2 



+9fta(3a) 3 - Srf(Ka) 3 - JJaSa (SrfSa + 3ta9fca) 



k 



*4 



(11) 



where 5(0) = or ±7r. Note however that, due to a rapid fall of the elasticity coefficient and 
the phase shift near the laboratory S-momentum of ~ 100 MeV/c, the use of the effective 
range approximation in Eq. is valid up to k* of ~ 50 MeV/c only. 

In Fig. [2, we show the iV£ correlation functions corresponding to isospin 3/2 (panel a) 
and 1/2 (panel b) as well as the singlet (S = 0) and triplet (5 = 1) contributions calculated 
according to Eq. ([6]) for the Gaussian radius ro = 3 fm. The enhancement and suppression 
at small k* is related with the positive and negative real parts of the scattering lengths, 
respectively. A wide suppression of the triplet contribution to the isospin-1/2 correlation 
function is due to large imaginary parts of the corresponding scattering length and effective 
radius. Though the effective range approximation is valid for this channel up to ~ 50 MeV/c, 
we do not expect a substantial change of the suppression form since at higher values of k* 
the correlation function already starts to approach unity. In any case, one may not rely on 
the effective range approximation in Eq. (Q and express the scattering amplitude directly 
through the elasticity coefficient and the phase shift according to Eq. ([H]). 
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FIG. 2: The NT, correlation functions corresponding to isospin 3/2 (panel a) and 1/2 (panel b) 
calculated for the Gaussian radius ro=3 fm assuming a uniform population of the spin states, i.e. 
po = 1/4 and p\ = 3/4. The singlet (S = 0) and triplet (S = 1) contributions are shown by the 
dashed-dotted and dashed curves, respectively. 

C. Two-channel approach 

The interaction of final state particles a and b can proceed not only through the elastic 
transition ab — > a& but also through inelastic reactions of the type cd — > ab, where c and d are 
also final state particles of the production process. The FSI effect on particle correlations is 
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known to be significant only for particles with a slow relative motion. Such particles continue 
to interact with each other after leaving the domain of particle production and their slow 
relative motion guarantees the possibility of the separation (factorization) of the amplitude 
of a slow FSI from the amplitude of a fast production process. For the relative motion of 
the particles involved in the FSI to be slow, the sums of the particle masses in the entrance 
and exit channels should be close to each other. Thus, in our case, one should account for 
the effect of inelastic transition nS + — > pE° in addition to the elastic transition pS° — > pL° . 
Instead of a single channel Schrodinger equation one should thus solve a two-channel one 
(the effect of the pA channel is taken into account in the complex effective single-channel 
.^/-functions in the isospin basis). In solving the standard scattering problem, one should 
take into account that the FSI problem corresponds to the inverse direction of time. As a 
result, one has to make the substitution k*(= k* = — k£) — > — k* and consider pS°(= 1) as 
the entrance channel and nE + (= 2) as the exit channel. Further, in single-channel equations 
©-(El), one has to substitute the amplitude f 3 , the i^-function, the low-energy scattering 
parameters a 3 , d 3 and the momentum k* by the corresponding symmetric 2x2 matrices 
f 3 , K s ', a 3 , d 3 and k: 

f s =[(K s y 1 -tky\ (K 3 )- 1 = (a 3 )- 1 + U s k* 2 . (12) 

The momentum matrix k is diagonal in the channel (particle) basis: kji = ki5ji] in accordance 
with the energy-momentum conservation in the transitions 1 — ■> i, k\ = k* a = kl = k* and 

1/2 



&2 — k c — k d 



*2 



2/U 2 f h m a + m b - m c - m d 



(13) 



where /xi = m a m b / (m a + m b ) and /i 2 = m c m d / (m c + m d ) are the reduced masses in the 
channels 1 = (a, b) and 2 = (c,d). Finally, the wave function ip 3 in Eq. ((Tj) should be 
generalized to the two-channel wave function vector %p s ' 11 describing the transitions 1 — ► i: 

tffcV) = exp(-ik«r-) + /f l( r)^ll, ^; (r .) = ^^SSmi, (14) 

r* 1 y yUj r* 

where r* = r* — r£ or r* — is the spatial separation of the particles in the exit channel. 

Since the particles in both channels are members of the same isospin multiplets, one can 
assume that they are produced with about the same probability. Therefore the correlation 
function will be simply a sum of the average squares of the wave functions ipj^l(r*) and 
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^-k*( r *) describing the respective elastic and inelastic transitions [9|. Similar to Eq. (1] 
one then has: 



ft(Pl,P2) 



(l^(r*)| 2 ) + (|^(r*)| 2 ) 

= 1 + 



5,21/ *\i2\ 



1 




2 


r 



2^(F) 
+ V^r l{Qo) 



To 



F 2 (Qr Q ) 



\- 1/^2 



/ 2 S i(^) 



r 



(15) 



where the analytical expression in Eq. (1151) corresponds to the Gaussian redistribution 
(|1J); since in our case the momentum k 2 > k± = k* is real, the contribution of the inelastic 
transition (the last term in Eq. (|15p ) merely coincides with the quadratic term in the 
contribution of the elastic transition after the substitution /f x — > (/i 2 / l^i) 1 ^ 2 fii- 

One should correct Eq. f[T5]) for the deviation of the spherical waves from the true scat- 
tered waves in the inner region of the short-range potential. The corresponding correction 
ATI is of comparable size to the effect of the second channel [9|. It is represented in a 
compact form in Eq. (125) of Ref. 8|, similar to the single-channel correction in Eq. (11 01) . 
In our case one has 



Aft(pi,p 2 ) = -{A^rlY^Ps + \flM 2 + 2K(/f 1 /*)d: 

5 



(16) 



where dfj = 2$ld(K s )^ / dk* 2 \ at k* = 0, d s coincides with the real part of the matrix of 
effective radii. 

Assuming that the isospin violation arises solely from the mass difference of the particles 
within a given isospin multiplet, one can express the elements of the matrices d s , d s , K s 
or (K 8 )^ 1 in the channel basis through the elements of the corresponding diagonal matrices 
in the representation of total isospin / (the products of the corresponding Clebsch-Gordan 
coefficients being 2/3, 1/3 and ±\/2/3). Particularly, 

2 



(^)n 1 



121 



5\-i 

1/2 



■;, _■- /3/2 (-^ 5 )l/2 



(17) 



Knowing the elements of the symmetric matrix [K ) , one can make the explicit inversion 



of the symmetric matrix (/ ) 1 given in Eq. (TT2"j) and get the required elements f$ of the 
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scattering amplitude matrix: 



DfZ = (f S rt 2 = (K S )-i 2 - ih Dfl = -(f s )- 2 \ = -(K s )- 2 \ 
Dfg = (f S rA = (K S rA-ik 1 

D = det^ 5 )- 1 = if) l: (.n, 2 - [(f S )~ii} 2 . (18) 

Note that at the momenta k\ = k* sufficiently larger than the momentum k 2 = 44.7 
MeV/c of the channel nS + at the threshold of the channel one can neglect the difference 
between the channel momenta and apply the relations 0171) directly to the elements of the 
amplitude matrix f s . 



D. Results 

The pS° correlation function as well as the singlet (5 = 0) and triplet (5 = 1) contribu- 
tions calculated for the Gaussian radius tq = 3 fm are shown in Fig. [31 As already mentioned 
in the discussion of Fig. [21 the enhancements at small k* are related with the positive real 
parts of the scattering lengths and a wide suppression of the triplet contribution is due to 
large imaginary parts of the isospin-1/2 scattering length and effective radius. 

In Fig. HI we compare the pA correlation function with the pS° and the residual pA^o 
ones calculated at the same conditions. Note that in our model the parent correlation 
functions are independent of the single particle spectra contrary to the residual correlations. 
To calculate the latter, we have used the thermal-like distribution ([2]) with To = 168 MeV/c. 
One may see that the pA^o residual correlation function is quite different from the pA one. In 
accordance with the discussion after Eq. (jHJ), the former is close to the parent pS° correlation 
function for k* > 70 Mev/c. In high energy heavy ion collisions the fraction of A's from 
S° decay is ~ 40%. If the corresponding residual correlation were neglected, the parent pA 
correlation effect would be underestimated thus leading to an overestimation of the source 
size. 

The residual py correlation functions resulting from the parent pS° correlation due to 
the S° — > A7 decay calculated for different Gaussian radii of the source are shown in 
Fig. [51 In accordance with the discussion after Eq. (jSJ), the parent correlation at small 
relative velocities is shifted to rather narrow fc*-region centered at ~ 70 MeV/c. Fig. [5] also 
demonstrates the sensitivity of the residual correlation effect to the source size. 
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FIG. 3: The pS° correlation function calculated for the Gaussian radius ro=3 fm assuming a 
uniform population of the spin states. The singlet (S = 0) and triplet (S = 1) contributions are 
shown by the dashed-dotted and dashed curves, respectively. 

It should be noted that there exists substantial uncertainty in the theoretical predictions 
for the low-energy scattering parameters in the isospin-1/2 iV£-channel. Thus the predic- 
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FIG. 4: The pS° (solid curve), the pA (dashed curve) and the residual pA^o (dashed-dotted curve) 
correlation functions calculated for the Gaussian radius rg=3 fm assuming a uniform population 
of the spin states and the thermal-like momentum distribution ([2]) with To = 168 MeV/c. The 
thin solid and thin dashed-dotted curves correspond to the p£° and the residual pA s o correlation 
functions calculated with the scattering parameters cl\j 2 = (2.54 + i0.26) fm, dy 2 = obtained 
from the pole position in the NSC89 model [l(J on the assumption of vanishing effective radius. 

tions of various Nijmegen potential models for the near-threshold pole position a^ 2 in this 
channel jlO] yield in the limit of zero effective radius, when a\, 2 = a\* 2 /\a\/ 2 \ 2 , similar triplet 
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FIG. 5: The residual j»7so correlation function resulting from the parent pT, correlation due to 
the S° — ► A7 decay calculated for the Gaussian source radius r$ = 2 fm (dashed curve), r$ = 3 
fm (solid curve) and tq = 5 fm (dashed-dotted curve). The thermal-like momentum distribution 
P| with To = 168 MeV/c is assumed for parent particles. The thin solid curve corresponds to 
ro = 3 fm and the scattering parameters ajy 2 = (2.54 + i0.26) fm, d\, 2 = obtained from the pole 
position in the NSC89 model [lo| on the assumption of vanishing effective radius. 



a 



scattering length a\, 2 as that deduced from Ref. |6j in the case of NSC97f and NF potentials 
while, they yield even opposite sign of the real part of this scattering length in the case of 
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earlier NSC89 and ND potentials ll|. To demonstrate the effect of possible uncertainty, we 
present in figures H] and besides the correlation functions corresponding to the potential 
model of Ref. 6j, also those obtained from the pole position in the NSC89 model assuming 
d\j 2 = [ill ]: a\, 2 = (2.54 -MO. 26) fm. One may conclude from these figures that the shape 
and the intensity of the pA s o and p^o residual correlations are sensitive to the pTP FSI and 
source size parameters thus providing a new possibility to learn about these parameters. 

The fraction of residual jry correlations arising from parent pTP correlations with respect 
to all other contributions into p'y system is not so large as for the pA system. The background 
arises mainly from photons from 7r° decay. Such photons can reduce the effect of our interest 
in pj system and make it invisible. The methods of the background suppression depend on 
experimental details and should be discussed separately. 



IV. CONCLUSION 



The two-photon and proton-photon residual correlations can serve as a new important 
source of information on the FSI and/or source size parameters as well as on the direct 
particle fractions. Particularly, a nontrivial femtoscopic irregularity in the proton-photon 
correlation function centered at k* pa 70 MeV/c is expected due to the P7so residual cor- 
relation. It is shown that the pA s o residual correlation recovers the negative parent pTP 
correlation function for k* > 70 Mev/c. The neglect of this negative residual correlation 
would lead to the underestimation of the parent pA correlation effect and to an overestima- 
tion of the source size. 
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